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Abstract: We analyze the effect of a finite volume on the thermodynamic potentials of a 
relativistic quantum field theory defined on a hypertorus at vanishing chemical potential. 
Using the symmetries of the Euclidean partition function, we interpret the thermodynamic 
observables as the expectation value of the energy-momentum tensor in the same theory 
living on a (3 x L 2 volume. In the case where the screening correlation lengths in the thermal 
system are finite, we obtain a closed formula for the leading finite volume effects in terms 
of the smallest screening mass. This formula can be used to estimate, and possibly correct 
for, the leading finite volume effects in lattice calculations of QCD thermodynamics. 
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1. Introduction 

The description of relativistic quantum systems at finite temperature plays a central role in 
cosmology, astrophysics, plasma physics and in the physics of heavy-ion collisions. In the 
latter context, the thermodynamics of Quantum Chromodynamics (QCD) is being studied 
intensively by lattice Monte-Carlo methods [1, 2, 3, 4, 5] and by analytic [6, 7] and semi- 
analytic methods [8]. The lattice calculations have to control both the discretization errors 
and the finite volume effects. References [9, 10, 11] address the question of how to reduce 
the former uncertainties. Here we address in some generality the finite-volume corrections 
to the energy density, entropy density and pressure calculated on a hypertorus of dimensions 
L 3 . At the same time, the method we follow provides a complementary point of view on 
the thermodynamics of the quantum field theory. This alternative interpretation might 
find some use in approximate analytic treatments such as the variational method [12, 13]. 

Finite size effects in gauge theories have been studied before at weak coupling [14, 
15]. In [16], an elegant calculation is presented that yields the finite-size effects for non- 
interacting gauge bosons. For non-Abelian gauge theories at extremely high temperatures, 
these finite- volume effects are indeed expected to be the leading ones in the regime 1/g ^> 
LT 3> 1 (L is the linear box size, T is the temperature and g is the gauge coupling). In 
that regime, the effects of electric and magnetic screening [17] are absent. However, at 
any finite temperature the asymptotic LT — ► 00 finite-size effects must be exponentially 



- 1 - 



suppressed by the finite, non-perturbative spatial correlation length. In other words the 
finite-size corrections to the thermodynamic potential is bound to be 0{e~ ca2TL ), where c 
is a number of order unity. At a few times the deconfining temperature T c , the screening 
masses are known to some extent from four-dimensional calculations, both in SU(3) gauge 
theory [18, 19] and in full QCD [20]. And at much higher temperatures, the dimensional 
reduction approach allows one to predict the temperature evolution of these correlation 
lengths [21]. Choosing L large compared to the longest of these correlation lengths should 
therefore ensure that finite-size effects are small. 

In this paper we first reinterpret the finite- volume effects by using symmetry properties 
of the Euclidean partition function and of the stress-energy tensor. We point out that 
there is a generic dynamical regime, where the leading finite size effects can be expressed 
completely in terms of the gap in the spatial screening spectrum and its derivative with 
respect to temperature. This is possible because T^ v plays the dual role of stress-energy and 
energy-momentum tensor: on the one hand its thermal expectation value gives the energy 
density and pressure, on the other hand its diagonal matrix elements on an individual 
state yield its energy and momentum. Since the screening gap can be calculated on the 
lattice relatively easily, the formula we derive allows one to estimate the finite-size effects 
in practice, and possibly to correct for them. 

In section 2 we describe the thermodynamic potentials as expectation values of ele- 
ments of the energy-momentum tensor and the dual interpretation of these matrix elements. 
In section 3 we exploit this interpretation further to estimate the leading finite-size effects 
on the thermodynamic potentials. Numerical applications to QCD are presented in section 
4, and we finish with some concluding remarks (section 5). 

2. Thermodynamic observables and their dual interpretation 

We consider a relativistic theory in four space-time dimensions without chemical potentials 
in Matsubara's Euclidean formalism. Euclidean expectation values are denoted by (.). At 
zero temperature and in infinite spatial volume, the system has a full SO (4) symmetry group 
corresponding to the Lorentz group in Minkovsky space. As a consequence of translation 
invariance, the theory possesses a conserved, symmetric energy-momentum tensor T^ v . We 
will sometimes consider separately its traceless part and its trace 9, using the notation 

= 6^ + -^95^, 9^ = 0. (2.1) 

The conserved charges measure energy and momentum, respectively, i.e. for a common 
eigenstate of these operators, we have 

J d 3 xf 00 (x)\V) =E\V), (2.2) 

J d 3 xf 0k (x)\V) =P k \V). (2.3) 

The partition function is Z = e~^™ in terms of the eigenvalues of J d 3 xToo, where 
(3 = 1/T is the inverse temperature. The pressure, energy density and entropy density are 
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obtained from Z according to 



p = T 



d\ogZ 

dv 



+cst, 



V 



d\ogZ 
dT 



+cst , 



s = 



v 



1 fd{T log Z) 
V { dT 



(2.4) 



v 



We will exclusively be considering volumes V = L 3 , with periodic boundary conditions in 
all four directions for bosons, and antiperiodic boundary conditions in all four directions for 
fermions. The energy density and pressure are defined up to an additive constant, which 
we choose such that both vanish at f3 = L. This is the standard choice in Monte-Carlo 
simulations. We have (e — 3p)((3,L) = (0)^x^3 — {9} L 4, (e + p)((3, L) = |(6>oo)/3xL 3 and, 
in the limit L — ► oo, s = (3(e +p). We remark that one can define different operators 
that play the role of energy-momentum tensor and lead to the same conserved charges [22]. 
But due to translation invariance, the Euclidean expectation values of the canonical and 
the Belinfante energy-momentum tensor are identical, since they differ only by a total 
derivative term [23]. 

In [24], using exact lattice QCD sum rules [25, 26, 27] we showed that if |^) is a state 
of definite energy E living in a periodic box L\ x L2 x L3, 



r 3 

<*| />*M*)|*> = \ i-^E L 

J L k=l 



d 
'dL k 



E. 



and one can similarly show that 



k=i 



d 
dLk 



E. 



(2.5) 
(2.6) 



l_4L 3 -^- + yL fe -^- 



k=i 



E. 



(2.7) 



The states are normalized such that (^l^) = l 1 . Strictly speaking, these equations hold 
when taking the difference between two states. Note however that for |^o) = the vacuum 
state in infinite volume, (^ol^ool^o) = (*o I ^33 1 ^0) = by Euclidean symmetry. Therefore 
the energy appearing on the right-hand-side of Eq. (2.5) and (2.7) can be thought of as 
the energy of the state |*) relative to the infinite- volume vacuum |*o)- For the operator 
8, due to its mixing with the unit operator one must always consider differences of matrix 
elements. Finally, we expect relations Eq. (2.5-2.7) to be true in other relativistic theories 
as well. 



2.1 Interchanging the coordinate axes 

Let us consider the expectation value of the operator #oo- Since dkk = —&oo, in a /3 x L 3 
box we have (8ij) = —^-(6*00)- We now want to reinterpret the axis labels 2 . The axis 3 
will play the role of Euclidean time (with extent L), while the short axis assumes the role 

1 In the infinite volume limit, one recovers from Eq. (2.5-2.7) the Minkovsky-space result (^IT^^I*) = 
P^Pv/M for covariantly normalized one-particle states, (^E , |^ r ) = (i5/M)L 3 . 

2 The idea of interchanging the coordinate axes in this way is of course not new, see for instance [28]. 
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of a spatial direction (with extent (3). In the expectation values below, we always indicate 
the dimensions of the lattice (in the order 0, 1, 2, 3). 

In this new system of coordinates, the operator #33 plays the role of the 00 component 
of the same tensor, 

(#Oo)/3xL 3 = -3 (#33)/3xL 3 = ~ 3 (6'oo)lx (L 2 /3) • ( 2 - 8 ) 

Next we apply formula (2.5) on the energy eigenstates of the (3 x L x L system. Note 
that for a homogeneous state, by which we mean E oc L1L2L3, this expression vanishes. 
However, we will apply this on the lowest-energy state of the (3 x L x L system. We write 
the energy levels of that system Eq, E±, E2 etc. ordered by increasing energy. We expect 
the energy per unit volume to have a finite limit when L — > 00, 



(2.9) 



The energy density eo(/3) is thus measured relative to the infinite-space vacuum. In this 
section we take the limit L — > 00 in Eq. (2.8) and assume that therefore the expectation 
value of a local operator is equal to its expectation value in the ground state of energy Eq. 
A sufficient condition for this is that there should be a spectral gap between E\ and Eq. 
By combining Eq. (2.5) and Eq. (2.8), we learn that the entropy density of the thermal 
system corresponds to 

s - -^\%)) / 3xoc3 - P ■ (2.10) 

Similarly, using Eq. (2.6) one easily finds that 

e-3p = 4e o ((3)+f3^0^ (2.11) 

By taking a linear combination of the last two equations, we also obtain the 'dual' inter- 
pretation of the pressure of the thermal system 

p = -eo(J3). (2.12) 

For instance, in a regime where the system behaves in a scale-invariant way, s = cT 3 and 
e — 3p = 0, the corresponding dual ground-state energy is given by 

W) = "^4- ( 2 - 13 ) 

We remark that finite-temperature phase transitions are mapped into quantum phase 
transitions in this interpretation [29]. The vacuum energy eo(J3) has a non-analyticity at 
a critical value of (3 equal to 1/T C . This non-analyticity is typically due to an avoided 
level-crossing. 
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3. Finite-volume effects on the thermodynamic potentials 



We can exploit the dual interpretation of the partition function further to study the finite- 
volume effects on the thermal system. Through a chain of relations, we successively relate 
the expectation value of the energy-momentum tensor on a [5 x L 3 lattice to the same 
expectation value on a j3 x oo 3 lattice. This allows us to arrive at a formula for the finite- 
volume correction to the thermodynamic potentials. Starting with the thermal expectation 
value of #oo> we successively write 

-g(#00)/3xL 3 = (#33)/3xL 3 = (#Oo)lx (L 2 x/3) = (#0o)oox (L 2 x/3) + ^1 ( 3 -l) 
= (#ll)Lx(ooxLx/3) + £l = ( 6l ll)oox(ooxLx/3) + K-i + )C 2 
= (#00)oox(ooxLx/3) + /Ci + /C 2 = (6 , 22)lx(ooxoox/3) + ICi + K, 2 

= (#22)oox(ooxoox/3) + /Ci + /C 2 + /C 3 = --(#0o)/3xoc.3 + /Ci + /C 2 + /C 3 . 

A spectral representation for the /Q is obtained in appendix A, for instance 

/Cl = -^TVTlT. {&n\fd 3 x M*)|*n> - <*o|/d 3 *M*)|*o>) e"^"^ , (3.2) 

where the |* n ) and E n are the eigenstates and energy levels of the (3 x L 2 system. We can 
use Eq. (3.1) to produce an expression for the finite- volume effects on the entropy density: 

4/3 4/3 

s = Y^°°hxoo 3 = y^ 00 ^ xL3 + ^(^i + ^2 + /c 3 ) . (3.3) 

Following the same steps as for the entropy density, we can obtain an expression for 
the leading finite-volume effects on the interaction measure. This case is slightly simpler, 
because the trace-anomaly operator is a Lorentz scalar: 

(#)/3xL 3 = (#)lx(L 2 x/3) = (#)oox(L 2 /3) + Jl (3-4) 

= (#)lx(coxLx/3) + Jl = (#)oox(ooxLx/3) +J1 + J2 

= (#)lx(ooxoox/3) +Jl+J2 = (0)oox(ooxoox/3) + Jl + J2 + Jz 

= (®}f3xoo 3 + Jl + J2 + J3 ■ 

A definition for the Ji based on the spectral representation is given in appendix B. We can 
apply the same reasoning to the L A system, sending the extent of each direction in turn to 
infinity. There are then four correction terms (2^) instead of three. Therefore we obtain 

e-3p= <0>/W-<0>oo4 = (9)/3xL3-(0)n-(Ji + J2 + J3) + (lo+li+l2+W- (3-5) 
3.1 The case of finite and discrete screening masses 

The general formulas (3.3) and (3.5) can be used together with the spectral definition of 
the TuiJiiK-i to predict the finite volume effects on the thermodynamic potentials. In the 
following, we make a qualitative assumption on the spectrum of the theory on a (5 x I? 
hypertorus with L>/?, which is in particular relevant to QCD at finite temperature. 
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We consider the case where the low-lying screening masses are discrete energy levels 
of the (3 x L 2 system. That is to say, 



Lm = L(E X - Eo) » 1 (3.6) 

and the next energy levels are simply that same excitation with non-zero momentum in 
the 'transverse' dimensions of size L, 



I Z7T 

w(kj_) = ym 2 + k2_, k ± = — (m,n 2 ), rii £ Z. (3.7) 



This lightest screening excitation can potentially have a u-iold degeneracy. The next screen- 
ing mass is assumed to be separated by a gap from the lowest one, [m^ — m)L 3> 1. In 
that situation, the /Q and Ji can be evaluated in a simple fashion, since they receive con- 
tributions only from one-'particle' states. We use particle in quotes because the lowest 
excitations have only two components of momentum; higher up in the screening spectrum 
one expects states with an additional energy of order 1/(3. 

We expect the scenario described above to apply in asymptotically free and conformal 
non-Abelian gauge theories. For every relativistic theory, the appropriate regime must be 
studied in order to correctly predict the leading finite- volume effects. As a counterexample 
to the above scenario, it is well-known that magnetic fields are not screened in an Abelian 
plasma. 

Since we are interested in the leading finite-volume effects, in the remainder of this 
section we write equations that hold up to terms of order max(e' 2mi , e~ m2L ). In appendix 
A, we calculate the corrections /Q = i 5 2,3 under these assumptions and find: 



ve -mL 



ve -mL 

£■2 = - 



3 mL 
2 + 2mL + -^m 2 L 2 —(3dpm((3) 

1 + mL + hm L + mL 



4 



(3.8) 
(3.9) 



2ir(3L 

Plugging these expressions into Eq. (3.3), we obtain our final formula 

4/3 mv e~ rnL 

s - y<0oo)/3xL3 = 2nL MP) ~ Wdpm{f3)] + ... (3.10) 

It shows that the knowledge of the longest spatial correlation length 1/m as a function of 
temperature T = 1/(3 allows one to compute the leading finite- volume corrections. 

The corrections 1^ and are computed in appendix B under the same assumptions 
formulated above. The zero-temperature volume corrections T M are assumed to be due to 
uq degenerate states of mass m^. We find to leading order 

— mL 

Ji = J 2 = Js = [m + (3dpm\ (3.11) 

J = Ji = J 2 = X 3 = Ki(m £) , (3.12) 
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where K\ is the modified Bessel function, rriQ is the mass gap of the theory on the hypertorus 
of size L 3 , and vq is its degeneracy. For instance, in isospin-symmetric QCD there would 
be vq = 3 pions. The final formula for the leading finite-volume effects on e — 3p follows, 

e ~ 3p = {6) pxL3 - <0) L 4 - [ml H + dpm] + Xi(m L) + . . . (3.13) 

Combining Eq. (3.13) and (3.10) with the thermodynamic identity Ts = e+p, we find that 
the pressure p is the thermodynamic quantity with the simplest finite- volume effect: 

P = ~\ ((Tkk)^ ~ (T kk ) Li ) + - gig ^(moL) + . . . (3.14) 

When the zero-temperature finite-volume corrections are negligible, the pressure computed 
in finite volume is lower than in the thermodynamic limit. Note that in Eq. (3.14) the 
pressure p(L) is assumed to be obtained directly from the expectation value of T k k- If p/T 4 
is obtained with the so-called 'integral method' (see for instance [30]), i.e. by integrating 
(e — 3p)/T 4 over temperature starting at T = 0, one should go back to Eq. (3.13) to 
compute the finite-volume effects. 



4. Applications 

We give two examples where we expect the formulas derived above to apply. 
4.1 Confined phase of SU(iV) gauge theory 

We first consider the pure SU(iV) gauge theory (see [31] for a review of its properties). 
Below the deconfining temperature T c , the center symmetry associated with the direction of 
length (3 is unbroken. Correspondingly, the expectation value of the Polyakov loop vanishes 
even in the infinite spatial-volume limit. However, for N = 2 and 3 the correlation length 
associated with the sector of non-zero winding number 3 becomes very long as the critical 
temperature is approached from below. In fact, in the case of SU(2) gauge theory, it even 
diverges with the 3d Ising exponent [32]. For SU(3), the correlation length becomes very 
long but remains finite. In [33], it was found that 

m(£)/T c = 0.53(4). (4.1) 

We can use the reinterpretation of the partition function to estimate the leading finite- 
volume effects on the pressure. At zero-temperature, the lightest state is the scalar glueball, 
so = 1. Given its large mass, Mq/T c ~ 5.3 [34, 35], it is not difficult to ensure that the 
corrections are negligible by making the box size L large enough. We therefore have 

p{T c ,L = c) p(T c ,L) 

YA = T 4 + 5 ' ( 4 - 2 ) 

3 For SU(3), the sector of winding number 2 is equivalent to the sector with winding number -1, which 
by charge conjugation has the same correlation length as the +1 sector. Therefore there are only two sector 
to discuss (winding and +1). 



-7- 



with 

m 2 e -{m/T c )-LT c 



T 2 2ttLT c 



= (0.0013, 0.00031) for LT C = (4, 6). (4.3) 



In fact, the value of p(T c )/T 4 is not known precisely, but is most likely on the order of 0.02, 
based on available numerical data [36], or on the pressure exerted by the known spectrum 
of glueballs [37, 38], assuming that they are non-interacting. Therefore the correction at 
LT C = 4 is not negligible if one aims at a precision of one percent on the pressure. 

In order to predict the finite volume correction to the entropy density, we need an 
estimate of the derivative of m with respect to (3 = 1/T. For a first idea of the order of 
magnitude involved, we can use the Nambu-Goto formula [39, 40, 41] for m([3) = a e f[(P)f3 
with 



2tt 1 



1 



(4-4) 



3 a/3 2 

where a is the string tension at T = 0, to estimate the derivative of the screening mass 
with respect to (3. The finite-volume effect is proportional to 

m(m - 3(3dpm) = -2(m 2 + ira) . (4.5) 

In particular, this quantity is negative, so the 'effective' entropy density computed in finite 
volume decreases towards the infinite-volume limit (the sign is opposite to the volume 
correction on the pressure). When m becomes small near T c , the magnitude of the finite 
volume effect on s is about ^j—. We can do a numerical application in the SU(3) case. 
For LT C = 4, using the value of m((3) given in Eq. (4.1) and T c /y/a « 0.64 [35], we get 
ae LT - A ~ 0.08. Since s/Tj? itself is about 0.2 [36, 42] (approaching from the confined phase), 
this is a large effect indeed. A box size of LT C = 9 is required to reduce this finite-size 
effect to one percent. This corresponds to a length L of about 6fm. 

It would be interesting to know in what range of quark masses this large finite- volume 
effect persists in full QCD, even though the center symmetry responsible for the existence 
of the light mode is badly broken in the presence of light quarks. Beyond checking that 
the leading correction is numerically small, it is also important that the exponent mL be 
large, as otherwise corrections that are formally higher order can be important. 



4.2 In the deconfined phase 

Let us consider the SU(3) gauge theory in the deconfined phase. Above T c , we know that 
the smallest screening mass corresponds to a state invariant under all symmetries of the 
theory in a [3 x L 2 box [18, 19] and its value is [43] 

f3m((3) = 2.62(16), 2.83(16), 2.88(10) (4.6) 

respectively at the temperatures 1.24T C , 1.65T C and 2.20T C . Due to these large values of 
the screening mass, the volume correction to the pressure appears to be negligible already 
for LT = 4, 8p/T 4 ps 8 • 10~ 6 . Recall that the Stefan-Boltzmann pressure is psb/T a = 
8vr 2 /45 ~ 1.75 for N = 3, N f = 0. 
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In the strongly coupled, large- N, J\f = 4 Super- Yang-Mills theory the screening masses 
have been calculated by AdS/CFT methods [44, 45]. They turn out to be significantly 
larger than in QCD, so that finite-volume effects would be even smaller for the same value 
of LT. 

At asymptotic temperatures in both QCD and SU(3) gauge theory, the smallest screen- 
ing mass corresponds to the Af + state of three-dimensional gauge theory, with a mass 
m/g 2 ~ 2.40 [21, 46], and g 2 = g 2 (T)T to leading order. When the coupling reaches the 
value it takes on the Z pole, a s « 0.11, this means that 5p/T 4 8 • 10 -7 for LT = 4. We 
conclude that the aspect ratio LT has to be increased only very slowly with temperature 
in order to accomodate the magnetic screening length 1/m ~ l/g 2 T. 

5. Concluding remarks 

We have derived a simple way to calculate the finite-volume effects affecting the energy 
density and pressure of a relativistic theory at zero chemical potential in terms of the 
spectrum of the same theory defined on a spatial hypercube with two large cycles and one 
of length (3 = l/T. When that spectrum is discrete, the leading finite- volume effects can 
be calculated completely in terms of the mass gap. It is almost obvious that the finite 
volume effects should be of order e~ mL , but we have shown that the prefactor is also 
entirely determined by the screening spectrum and its temperature dependence. This is 
because the diagonal matrix elements of the energy-momentum tensor are themselves given 
in terms of that spectrum (see Eq. (2.5-2.7)). 

It is hoped that Eq. (3.10) and (3.13) will be useful in controlling the finite- volume 
effects in lattice QCD thermodynamics calculations. If the screening mass gap is known and 
a finite-volume study shows that the finite-size effects are well described by the formula, 
one can use it to correct for these effects. If there are several screening masses below the 
threshold of 2m(/3), those states will contribute terms to the finite-size effects similar to 
the lightest one. To include the effects of screening states above 2m((3), one presumably 
needs to know the scattering length of these 'particles'. In QCD at low temperatures, 
explicit calculations using chiral perturbation theory are then likely to be predictive, since 
information on the scattering lengths of pions is available. 

It is clear that the method followed here is not specific to four dimensions. It also 
applies for instance to three-dimensional gauge theories. The main difference is that the 
transverse momentum of the lightest screening state only has one component, so that 
momentum integrals as in Eq. (A. 6) become one-dimensional. 

In SU(iV) gauge theories, it is interesting to note the dependence of the finite vol- 
ume effects of energy density and pressure on the number of colors N. We showed that 
the asymptotic finite volume effect is driven by a unique color singlet state and that its 
contribution is therefore 0(N°). In the deconfined phase, the thermodynamic potentials 
are 0(iV 2 ), and the relative size of finite-volume effects is thus 1/N 2 suppressed. This 
conclusion remains qualitatively valid in the presence of quarks, since their main effect is 
to add a contribution of order NNf to the thermodynamic potentials. 
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In the confined phase on the other hand, the thermodynamic potentials are 0(N°), 
so there is no parametric suppression of the volume effect there. Since on the lattice the 
entropy density is simply computed as the difference between the lxl electric and magnetic 
Wilson loops ('plaquettes'), this might seem to contradict the statement that finite volume 
effects on expectation values of Wilson loops vanish in the large- N limit as long as the 
center symmetries remain intact. However the statement of volume- independence only 
applies to the 0(iV 2 ) contribution to the plaquette (see section 2.3 of [47] for a clear 
discussion); the latter is divergent in the continuum limit and cancels in the difference of 
electric and magnetic plaquettes. Therefore there is no contradiction between our results 
and the large- N volume-independence arguments. 
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A. Calculation of /Ci i2 ,3 

In this appendix, we compute the quantities /Q that are the finite-time extent corrections 
to certain expectation values of 9^ u . Let us start with JC±, it is the difference between the 
expectation value of (#oo) on a L x (L 2 x (3) lattice and on an oo x (L 2 x (3) lattice. (In 
this appendix, we drop the ~ on the energies and states of the (3 x L 2 system, and we set the 
degeneracy v of the energy level E\ to one, since the more general result is simply obtained 
by multiplying /Ci by v.) 

^ 2 ^00)lx(L 2 x/3) = Z{L \ L ^j E e ~ £ " L ^l / d 3 X§oo(x)\* n ) (A.l) 

= <*o| J d 3 x6 00 (x)\* ) 

+ (<*»l/^floo(*)l*n> " <*o| P 3 xM*)|*o>) e- (Mo)L . 

n>l 

We now observe that 

oox(L 2 x/3) • 

(A.2) 

We therefore identify JC\ as (see Eq. (3.1)) 



£i = J d 3 x§ 00 (x)\* n ) - (* | j d 3 x^oo(x)|*o)) e 



-(E n -E )L 

(A.3) 

We now assume that the lowest-lying excited states are 'one-particle' excitations with 
arbitrary momentum in the two directions of length L. Let us therefore call m = E\ — Eq 
the energy of the first excited state, since it is at rest. By rotation symmetry among the 
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three dimensions that are not of length (3, the dispersion relation must be relativistic and 
we define 

w(k) = Vm 2 + k 2 . (A.4) 
Here k is a two-component vector. The partition function is for instance given by 

Ze E L = £ e -(E n -E )L = 1 + ^2 e~" (k)L + 0(e~ 2mL ) . (A.5) 

n k 

but to the order we are working at, we can use Ze E ° L = 1. We are thus lead to the 
expression 

Kl = Jf 2 E e "" (k)L ((*i(k)|/d 3 x^oo(x)|*i(k)) - ^o\fd 3 x§ 00 (x)\^o))+... (A.6) 
Using Eq. 2.5, we obtain 

<*I^k)| fd 3 *M*)l*i(k)> - (* \fd 3 x8 00 (x)\* ) = w (k) - i __( m /3) . 
Using the Poisson summation formula and performing the integral, one finds 



/Ci = 



1 



E 



- \mdp{m0) 



L 4 m 2 - y 2 (1 + mVy 2 + L 2 ) + L 2 (2 + m 2 y 2 + 2m yV + L 2 ) 



(A.7) 



(y 2 n + l 2 ) 2 



Here y n = Ln and n £ Z 2 . It is now obvious that the terms with n / are subleading 
and can be dropped, hence 



-mL 



/Ci 



2vr/3L 3 
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2 + 2mL + m 2 L 2 - ±mL 2 — (m/3) 
4 op 



+ o( e - 2mL ) 



(A.8) 



Since we are assuming that mL 3> 1, we can always replace the momentum sum for a 
direction of size L by an integral, j^^2 k —> J f^- The corrections to this are suppressed 
by e~ mL , as just shown. It is therefore clear from their definitions that /C2 = K.3 to this 
accuracy. These expressions are calculated in the same fashion as K.±, using this time 
Eq. (2.7). The contribution of the |^i(k)) states reads 



m 



d 2 k 



,-w(k)L 



4/3 J (2tt) 5 



2c<j(k) (3d pm — m 



m 



W (k) 



+ 



(A.9) 



which leads to Eq. (3.9). 



B. Calculation of ^1,2,3 

The derivation follows closely that of appendix A. The spectral representation for J\ is 
1 



Ji 



f3L 2 Ze E » L 



((*n\fd 3 x§(x)\y n )} - (* |/d 3 x^>)|^o)) e-^~ E ^ L . (B.l) 



n>l 
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We use Eq. (2.6) to reach the expression 



Jl = j{m + Pa0m) j——— (B.2) 

Performing the momentum integral leads to Eq. (3.11). It is clear that the J. L only differ 
by the size of the spatial dimensions transverse to the dimension of size (5 (they are either 
of size L or infinite). Since we assume mL S> 1, this difference is sub leading, as can easily 
be seen by using the Poisson summation formula. 

As for 1 , the same definition as Eq. (B.l) holds, except that the states live on an 
L x L x L hypercube. Therefore we get (for a degeneracy uq of 1 and setting oo(k) = 
a/Ic 2 + m.Q, where mo is the mass gap) 

For bosonic degrees of freedom, the same calculation applies to the Ik, neglecting terms 
of order e ~ 2m o L m Therefore all Z u are equal. Indeed all directions are truly symmetric if 
the boundary conditions are periodic. For fermions, strictly speaking the boundary 
condition has to be antiperiodic in all directions for the same to apply. Indeed it is for 
antiperiodic boundary conditions that the path integral computes the trace over states. 
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